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ABSTJ^ 

A severe mismatch hetween generation and demand 
caused hy a sudden change of generation or load can cause the 
power system frequency to deviate appreciably from its normal 
value. The system frequency , could fall or increase ^ dependent 
upon whether the generation or the load is lost. Of course, 
the most critical case is the loss of one of the larger 
generating units of the system. This will load to under 
frequency :;onditions which cannot be tolerated for a long 
time even if the drop is about 3/» from the normal value. 

If the spinning reserve is also available, the frequency 
may recover within an acceptable time, otherwise load 
shedding may have to be done. Obviously, in order to evolve 
a proper operating strategy, we need to obtain the power- 
frequency tranvSient upon loss of generation. A common method 
is by extended transient stability studies. However, such 
simulation can be computationally very expensive because of 
small time step involved. In 1971, Stanton proposed 
'dynamic energy balance concept’ for sinKilation of power- 
frequency transients. In this method by avoiding synchroni- 
zing oscillations and other small time constants involved 
with subtransient phenomenon,- he claimed a ten to one speed 
advantage over the conventional transient stability method. 



Xlll 

'fhe aim of this thesis is to review the Stanton's 
scheme in the context of a practical problem faced hy a 
typical electricity hoard of our country and to develop a 
model of our own. fhe test problem chosen is taken from 
U.P. State Electricity Board planning model network of 
year 19S4-S5. 



OHAPTER 1 


IITTRODUOTIOR 


1 .1 General Discussion 

In power system analysis, the study of 'Power- 
frequency transients' has gained importance in order to 
analyse hehaviour of various power system components and 
also the system's "behaviour on the whole under abnormal 
conditions of frequency excursions in the system. It is 
a well known fact that the frequency of the system is 
very much dependent on the power balance between genera- 
tion and demand. Loss of a large block of generation or 
load will cause a severe mismatch between generation and 
demand leading to dangerously low or high frequency 
conditions in the system in few seconds after the occiarrence 
of the disturbance. Under these abnormal frequency condi- 
tions, sane of the electrical components like synchronous 
generators, large capacity synchronous motors etc. which are 
sensitive ,to frequency changes may loose synchronism causing 

f 

disastrous conditions in the systan. Mechanical vibrations 
of considerable amplitudes will be induced in the turbine 
shafts and blades due to the coincidence of their natural 
frequency modes with the system low frequencies following 
a loss of large block of generation causing permanent 
mechanical deformations in the shafts and blades. The modem 
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turbines must be run at speeds which are very close to 
synchronous to avoid permanent damage. Por example, 
specifications of 220 m machines which are being installed 
in our country now state that the turbines may be damaged 
permanently even if the speed corresponds to 48 Hz for few 
minutes. 

Following a loss of large block of generation, the 
redistribution of powers in the network may overload some 
of the lines already operating at their normal capacity 
causing line trip out conditions posing a threat to the 
stability of the system. To avoid such unpleasant 
situations, the system frequency should be brought back 
to its normal operating value in minimum amount of time, 
large frequency dips can be minimized by fast acting 
gpeed-governor mechanisms, by providing spinning reserve 
in the system’ and by shedding off the loads using fast 
acting frequency sensitive relays, fhe most common 
practice of restoring system's frequency under these 
transient conditions is by load shedding. Knowledge about 
capability of spinning reserve, behaviour of loads at 
abnormal frequency and voltage conditions will help in 
developing the most effective load shedding strategies. 

Since the disturbance location is most random 
in nature, exhaustive off-line power-frequency transient 
studies need to be carried out to develop efficient and 
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effective load shedding strategies so that minimum number 
of power consumers are affected. This can involve large 
amounts of computation effort resulting in excessive costs. 

So it is highly desirable to have methods which simulate 
power-frequency transients in an economical way, 

1 . 2 Previous Work 

The power-frequency transients resulting from 

severe mismatches between generation and demand can be 

simulated by extended transient stability methods. The 

inclusion of generator subtransient phenomenon, damper 

bar windings, small excitation time constants set a very 

low upper limit on time step of integration. However, this 

can be increased by appropriate choice of dynamic models. 

Honoe the ultimate limit on step size in transient stability 

studies is the ’frequency of synchronizing oscillations’. 

PE RtOD 

However, the li wratl#*! of these oscillations is also very 
small, a very low step size for integration cannot be 
avoided. Obviously this results in excessive computation 
costs for studying power-frequency transients. Even though 
extended transient stability studies do simulate power- 
frequency transients, it is primarily iieant for s-imulation 
of power angle transients which is not our aim. 

Rooently a new technique of simulation of power- 
frequency transients has been developed by Stanton [1] using 
’lynamio Energy Balance’ concepts. The development of dynamic 
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energy balance studios is motivated by the need for more 
economical dynamic simulation of power -frequency transients. 
The key steps involved in simulation procedure are elimination 
of synchronizing oscillations, small time constants used in 
representation of subsystems like excitation, speed governor, 
turbine and generator models. 


IXiring transient conditions, the entire system is 
assumed to have an equivalent speed which will be 

equal to the inertia weighted average of individual generator 
rotor speeds. This equivalent speed of the system is defined 
as the velocity of ’system centre of inertia’ . The centre of 
inertia is the rotational analog of centre of mass and is 
referred as centre of angle ’9^'* centre of angle 0^ is 
defined as 


9. 


_1_ 

M, 


n 

I 


T 


M. 9^ 
D 0 


( 1 . 1 ) 


where M., 


and Mgi 


moment of inertia of D'th machine, 

angle of the jth rotor measured with respect to 

a synchronously rotating frame of reference 

n 

t is the total moment of inertia of the 

3-1 ^ 


system. 

n «. number of machines present in the system. 

The dynamic equations of motion of system centre of inertia 
are described by the following equations. 
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« “o 1 , 

and 

^ ®o 

-ST = (“o - “s) (1-5) 

■where =? Total mechanical power input to the system, 

^eT “ ^otal electrical power output of the system 
and Wg = Synchronous speed of the system. 

The frequency 'f’ obtained by the solution of equation (1.2) 
will be the inertia weighted average of frequencies of 
individual rotors present in the system. The above equations 
form the basis of simulation of power -frequency transients 
by dynamic energy balance concepts. 

The synchronizing oscillations are eliminated from 
the simulation method by assuming that the oscillatory 
power transfers from one rotor to another during these 
oscillations are zero. These oscillatory power transfers 
do not significantly effect the total electrical power 
output Pgg, of the system. However, thqy may increase losses 
Slightly. Assumption of Pg^ constant during these oscilla- 
tions gives a scope for elimination of synchronizing oscilla- 
tions from the model. Elimination of these oscillations allows 
us to choose large time steps for integration and this 
facilitates to study the slower dynamics of the power System. 
It has been claimed by Stanton that the simulation of 
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p ow er -fr equ ency transients by this method requires one tenth 
the computational effort of a transient stability study yet 
giving similar results within the practical limits. 

1 .3 Details, of Problem under Oonsideration 

It appears from the literature that only Stanton 
has made a serious effort to develop a method for the study 
of power-frequency transients. The aim of this thesis is to 
review this method in the context of a practical problem 
faced by a typical electricity board of our country and 
develop a suitable model of our own. 

The problem chosen is taken from U.P. State 
Electricity Board planning model network for the year 
1984-85. ^ single line diagram of the planned network is 
given in K.g. 1.1.* It consists of 20 generators, 96 buses 
and 136 lines. The number of buses is not too large as only 
high voltage network is considered. Out of these 20 genera- 
ting stations, there are 11 thermal and 9 hydro stations. 

The total MYA capacity of the system is 7576 out of which 
1790 MYA is hydro and the rest, i.e. 5786 MYA is thermal. 

The objective is to study the power -frequency transient if 
52 s MYA thermal generator at Anpara is tripped. This 
particular generator, prior to the trip out is assumed to 
be delivering 400 MW power. The total spinning reserve ■ 
available at the time of th e trip out is assumed to be 
560 MYA. 

*The bus data prior to disturbance, line, transformer, 
shunt load and voltage controlled bus data is given 
in Appendix I. 
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The chapterwise description of the thesis is 
given helow. 

The second chapter deals with the speed governor 
and turbine models. 

The third chapter covers the simulation of power- 
frequency transients treating electrical loading as 
c on stan t . 

The fourth chapter deals- with the generator and 
excitation models. 

The fifth chapter deals with the simulation of 
power-frequency transients considering the 
electrical loading. 

Overall conclusions and suggestions for future 
work are projected in the last chapter. 
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CHAPTER 2 

TUHBIHE ARE SPEED GOVERNOR MODELS 


2 . 1 General Discussion 

Before we can study the power-frequency transient 
problem of the system described in Sec. 1.3, we need to get 
acquainted with the models for hydro and steam turbines 
and speed governing mechanisms. Therefore, they are 
discussed first as given below from Ref. 12], 


2 . 2 Speed Governor and Thermal Turbine Model 

The block diagram model representation of speed 
governor and thermal turbine are given in Pig. 2.1, In the 
speed governor block, the meaning of various quantities 
is as follows: 

^ref “ Heference speed, rad/ sec. 
oj = Actual speed, rad/sec. 

R = Speed regulation constant of speed governor 

f^ = Normal operating frequency, Hertz 

T = Speed governor time constant, seconds 
s 

= Lower limit on rate of change of power imposed by 
Jontrol valve rate limits 

^ax ~ limit on rate of change of power imposed by 

control valve rate limits 
g^^^ = Lower limit imposed by valve travel 

IX 


g ^^ ^ = Upper limit imposed by valve travel 
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g = Valve position, p.u. 

Similarly the meaning of various symbols used for the turbine 
block are 


= Time constants of high pressure, intermediate 
pressure, low pressure stages respectively in 
seconds 

Php> ^jp> P]^p=5'2^ac'tional power coefficients of high pressure, 
intermediate pressure and low pressure stages 
respectively 




'base 


base 


-m 


= 1 

= Turbine MVA rating 
= System base MVA 
= Mechanical power, p.u. 


The analog diagram representation of Pig. 2.1 is 
shown in Pig. 2. 2. The analog representation helps in 
writing the appropriate dynamic equations for digital 
computer solution as given below. Por calculation of 
initial conditions, the following equations are used. 


S = ^^m*^base^/^base*^^HP'^^IP^^IP^ 

(2.1) 

ref =0,0 . (g.R.2irfe) 

(2.2) 

^4= g 

(2.5) 

X5= g 

(2.4) 

Xg= g 

(2.5) 

X^= g 

(2.6) 



During step ty step integration procedure, the folio-wing 


equations a^e, used 

S - ^4 * ^in — ^4 — ^ax 

" ^“ref - )/(^^-2»f ) (2.S) 

^2 = - g)/Tg (2.9) 

^3 ~ ^2 ’ ^in — -^2 — ^ax (2.10) 

^4 = X3 (2.11) 


Equations (2.7) to (2.11) discrihe the performance of the 
speed governor. 

The solution of differential equation (2.11) gives 
the position of the effective controlled governor valve 
symbolised as g. This becomes the input to tho -turbine. 


The turbine equations are 



— (g ~ HP 

(2.12) 


ig = (X 5 - Xg)/Tjp 

(2.13) 


= (Xg - 

(2.14) 

amd 


(2.15) 


In the above equations, the differential equations are 
represented as X = dX/dt. The differential equations (2.11) 
to ( 2 , 14 ) are solved at each integration time step by making 
use of equations (2.7) to (2.10) to get the mechanical power 
output of the -turbine 















ilNE . . . 

SPEED QOve.RNOK 




r, ?f<) 


'4,V, 




''. ' " it' ' 1' 5 

r/^ 




M 


v.n^V 


v;; ^fr\f 
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2. 5 Speed Groverrcor and Hydro Turtine Model 

The hlock diagram representation of the hydro speed 
governor and turbine are given in Pig. 2.3. In the speed _ 
governor block, the various symbols have the following 
meaning. 

“ref’ ^s’ ^min’ ^ax’ ^in’ ^ax> S 
have similar meaning as given in Sec. 2 . 2 , 

Tj^ = Dash pot time constant, seconds 
o = Permanent speed droop coefficient 

6 = Transient speed droop coefficient 

Similarly the meaning of various symbols used in the 
turbine block are 

= Water starting time, seconds 

P, = Turbine liVA rating 
base 

S, = System base MVA chosen, 
base 

The analog diagram representation of Ptg.2.3 is 
shown in Pig. 2.4 from which the following equations are 
written. The initial conditions are calculated by the 


following 

equations . 


X/ 

= P„ ^base 

(2.16) 

4 

base 



- X ^ 

(2.17) 

Xy . 

= 1.5 X 4 

(2.1S) 

“ref 


(2. IQ) 
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Bering step by step integration, the following equations 
are used 


S ^4 » ^in — ^^4 — ^ax 

^1 " ^“ref " ‘*')/2Tif^ 

= g. 6 - 


Ig = (IC, - Xg - g .a )/5)g 

^5 = ^2 » ^in^ ^2 .< ^ax 

Xg = 2.0(x^-g) 


X 


X, 


4 - -3 


R 


X7 = (g - x«)A' 


¥ 


"" ^8 S. 


base 


base 


( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 

(2.24) 

(2.25) 

(2.26) 

(2.27) 

(2.28) 
(2.29) 


In the above equations the differential equations (2.26) to 
(2.28) are solved at each integration time step by making 
use of equations (2.20) to (2.25) to get the meohanical 


power output of the turbine, P 


m 


2 . 4 Turbine and Speed Governor I>ata for the UPSEB System 

The system has 20 machines out of which the turbine 
data is available for only 16 machines. The remaining 4 are 
treated as constant mechanical power output turbines through- 
out the simulation period. 'it- 
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The data for thermal turbines and speed governors is 
given in Table 2.1 and the data for hydro turbines and 
speed governors is given in Table 2.2. 



Table 2.1: Thermal Turbine and Speed Governor Data, 
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Table 2.2: Hydro Turbine and Speed Governor Data. 
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OHAPTER 3 

IMPIEMERTATIOR TREATING ELEDTRIOAI LOADING AS gONSTANT 


3.1 General Elscussion 

To get a ’feel* of the prohlem posed in Sec. 1.3, 
it was decided that at first we would not get involved with 
the electrical network and generator characteristics. Thus 
for the time being we shall assume that after the disturbance, 
the electrical loading is not changed and therefore the total 
load requirement plus the network losses stay constant. The 
moment the disturbance occurs, obviously the total output 
of the generation is less than the load requirement. This 
will result in lower speed than the synchronous value. 

This would inturn activate the speed' governor mechanism 
to increase the mechanical power output of the turbines 
whenever spinning reserve is available. Since in the 
problem under consideration, adequate spinning reserve 
of 560 MVA compared t-' the generation drop of 400 MW is 
available, eventually the system generation will match the 
system load plus the network losses. This problem was 
analysed in four parts as follows. 

Part -(a) ; It was assumed that the entire generation is 
thermal and it would be represented by one equivalent 


thermal turbine. 
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Part - (b) : It was assumed that the entire generation is 
hydro and it could he represented hy one equivalent hydro 
turbine. 

Part-(c ) : The ^stem is represented hy both the equivalent 

thermal and hydro turbines. All the hydro stations are 
represented by a single equivalent hydro turbine and all 
the thermal stations by a single equivalent thermal turbine. 

Part - (d) : In this part of study, each power station is 
represented by an equivalent hydro or thermal turbine as 
the case may be. Identity of all the power stations in the 
syston is retained. 

3. 2 Implementation 
Part - (a) : 

Ihe M7A rating of the equivalent thermal turbine is 
assumed to be equal to the sum of MVA ratings of individual 
power stations and is equal to 7576.25. It "was observed 
from the data given in Table 2.1 that all the speed govemers 
of the machines have similar characteristics, and also for 
the most of the turbines, the characteristics of high 
pressure, intermediate pressure and low pressure stages 
are similar. Prom these observations, the following data 
has been chosen for the equivalent thermal turbine. 

P^ = 7565.25 MYA 1^ = 0.2 seconds 

ba se , 

= 4.0 seconds = 1000.0 seconds 
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11 

0.62 

Ejp 

= 0.33 

i! 

0.0 

R 

= 0.04 

T = 

S 

0.425 seconds 

• 

^in 

= -1.13 

• 

^ax ” 

0.59 

^in 

= 0.0 and = 0.3 


Prom the initial load flow conditions prior to the 
disturhance, the total electrical power output of the system 
is calculated hy 


P 


el 


n 

I ^ • 

D=1 


(3.1) 


and is equal to 5556.37 MW 


P . = Electrical power output of the jth machine 


n 


Wumher of machines = 20 


Similarly the total mechanical power input to the system is 
computed hy 


n 

p = y P - 

ml m3 

where 

The system equations of motion are given hy 
d w ^ 

It " Mgn'^^mT ^eT^ 

where and are defined in Sec.1 .2 and 

^mT 


(3.2) 


(3.3) 




(3.4) 
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The equivalent thermal turbine dynamical equations described 
in Sec. 2.2 are "written symbolically in the form of equation 
(3.4). During each time step, the equations (3.5) and (3.4) 
are solved together by using Runge-Kutta Fourth Order method 
to obtain the frequency of the system given by the equation 

f = (3.5) 

The frequency is plotted against time and is shown in 
J'ig.3.1. 

Part - (b) : 

The following test data has been chosen for the 
equivalent hydro turbine on the basis of the similar 
observations discussed in part (a). 


p 

base 

= 7576.25 MVA 

T 

1.63 

seconds 


= 0.015 seconds 

T = 

1r 

12.0 

seconds 

a 

= 0.04 

5 

0.36 


^in 

= -0.2 

• 

^ax 

0.1 


^in 

= 0.0 

^ = 
‘^max 

0.9 



The simulation procedure is similar as given in Sec. 3.2, 

Part (b). The dynamical equations of hydro turbine and speed 
governor are described in Sec. 2.3. 

The frequency versus time graph is sho"wn in Pig. 3. 2. 
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Part - (o) ; 

The MYA rating of the equivalent thermal turbine is 
assumed to he equal to the sum of MYA ratings of the indivi- 
dual thermal power stations and is equal to 5786»0. Similarly 
the MYA rating of the equivalent hydro turbine is equal to 
the sum of MYA ratings of the individual hydro power stations.* 
The total load on the systan is 5556.87 MW. The same data 
given in Part-(a) and Part-(b) is used for the equivalent 
thermal and hydro turbines respectively. Prom the initial 
load flow data prior to the disturbance, total electrical 
power outputs of thermal and hydrogenerators are computed 
and are symbolised as • ''^ow the total 

electrical power output of the system is 

^eT " ^eT(Th) * ^eT(H) , 

Similarly the total mechanical power input to the system 
is given by 


^mT " ^mT(Th) * ^mT(H) 


(5.7) 


The equations of motion of the system are 


d“ 

c 

dt 


1_(P _ p ) 

mT eT^ 


^mT(Th) ^^^mT(Th)* 


^mT(H) = ^(^mT(H)’ 


( 3 . 3 ) 

(3.9) 

( 3 . 10 ) 


where equations (3.9) and (3.10), are symbolic representations 
of the dynamic equations of thermal and hydro turbines given 


* and is equal to 1790 MYA. 
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in Sec. 2.2 and Sec. 2. 3 respectively. At each time step 
of integration , equations (3.S) to (3.10) are solved 
together using B-unge-Kutta Pourth Order integration method. 
The graph of frequency versus time is shown in Pig. 3. 3. 

Part - (d) ; 

Por this case study , the turbine data used is 
given in Tables 2.1 and 2.2. Since all the turbines are 
represented by their respective speed governor controls, 
the dynamic equations of speed governor and turbines are 
written in the form 

“o^ ^here d = 1,2,...,n (3.11) 

equations of motion of the system are solved at each time 
step by using Runge-Kutta Pourth order integration method. 
The frequency is plotted against time and is also shown in 
Pig. 3.3. 

3 .3 Piscussion and Qonclusions 

It should be apparent that the results obtained 
in Part-(a) and Part-(b) cannot be close to the actual 
results. Representation of the system by one equivalent 
thermal or hydro machine was carried out in order to 
examine the relative performances of a purely thermal and 
a purely hydro systm. It is observed from Pigs. 3.1 and 
3.2 that a thermal system can pick up the load much faster 
compared to the hydro system. As seen from Pig.3.1, for a 
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thermal system, the frequency dropped only hy about 0.62 Hz. 
The minimum frequency condition occurred in less than 4 seconds 
and reached the steady state value in about 16 seconds. On 
the other hand, for the hydro systm, the frequency had gone 
down to as low as 47 Hz in less than 4 seconds and the 
oscillatory conditions continued to prevail as long as the 
duration of the study which is 40 seconds. This is because 
of high moment of inertia of the hydro turbines and also 
large time constants involved in the system dynamics. 

If we examine fig. 3. 3, it is seen that the behaviour 
of frequency of the system is in between figs. 3.1 and 3.2 
with no oscillatory behaviour. It was also observed that 
the representation of the entire system by one equivalent 
thermal and one equivalent hydro machine gave results very 
close to the one where the identity of all the power 
stations were retained (Part - d). This is due to the 
fact that most of the thermal as well as hydro machines 
have similar characteristics. 
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OHilPTER 4 

GmERATOR AWP EXCITATION SYSTEM MODELS 


4.1 Generator Representation 


A simplified generator representation is used Ly 
eliminating damper bar windings and other subtransient 
effects l3 1. Because of elimination of synchronising 
oscillations, which is di scribed’ in Sec. 1.2, relativ'e 
movement between the rotor and airgap flux wave is not 
rapid enough nor of sufficient magnitude to induce currents 
in the damper bar windings. Hence the damper bar winding is 


not represented in the model. The changes in field flux 


linkages are taken into account.With conventional d-q 
representation of the machine, only field winding is 
assumed to be on the d-axis. The quadrature axis is chosen 
to lead the direct axis by 90 electrical degrees. The 
position of the quadrature axis is determined by calcula- 
ting a fictitious voltage located on this axis. This is 
a voltage back of quadrature axis synchronous reactance 

X and is determined by the equation (4.1) 

Q. 



= + 


(r^ ^ 


(4.1) 


where E^^ 


Voltage behind quadrature axis synchronous 


reactance. 


= Terminal voltage of the generator. 
I^ = Terminal current of the generator. 
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= Stator resistance 

snd = Quadrature-axis synchronous reactance. 

The equivalent circuit of the generator and the phasor 
diagrams for steady state conditions are shown in fig. 4. "I. 


The voltage induced hy the field flux will he in 
phase with the quadrature axis and this voltage can he 
determined hy the equation 


where and are direct and quadrature axis components 
of terminal current respectively. The quadrature axis 
component of voltage haoLc of transient reactance is 
given hy 


E' 

q 


\ - 3(V^a>^d 


( 4 . 3 ) 


where E^ 


X’. = 


Voltage proportional to field flux linkages 
resulting from the combined effect of field 
and aima ture currents. 

Direct axis transient reactance and 
Direct axis component of terminal (current 


Since field flux linkages do not change instantaneously, 
E’ also does not change instantaneously. The rate of 

q ' , ^ 

change of E’ along the quadrature axis is dependent on the 

q 

field voltage controlled hy the regulator and exciter, the 
voltage proportional to field current and- direct axis 
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"transient open cirouit time constant and is given by the 
expression 

d. E ^ 

■^do 

■where Ej, = Eield voltage acting along the q'uadrature axis 

Ej = Voltage proportional to field current and 

= Eirect-axis transiait open circuit time constant, 
seconds. 

The detailed phasor diagram describing E , and E’ is 
give‘s in Eig. 4.1. 

El = E- . (X^ - X|)l3 (4.5) 

and Eq^ = E^ + (4.6) 

The above equations (4.4) to (4.6) represent the generator 
behaviour under transient state. The transient state vector 
diagram is shown in Eig. 4. 2. 

4 . 2 Ex:iitation System Model 

The excitation systan is represented by IEEE type-1 
rotating exciter model [4]. The block diagram is shown in 
Eig. 4.3. Various symbols used in the block diagram carry 
the following meanings. 

Vj^g^ = Regulator reference voltage setting, p.u. 

Vg, = (Generator terminal voltage, p.u. 




I 



FIG-4'2 VECTOR DIAGRAM OF SYNCHRONOUS GENERATOR 
IN THE TRANSIENT STATE 
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= Regulator gain 

^A ~ Regulator amplifier time constant, seconds 
Vr = Maximum value of regulator output voltage 

^Rmin = Minimum value of regulator output voltage 

% = Exciter constant related to self excited field 

Tj, = Exciter time constant, seconds 
Kj, = Regulator stabilizing circuit gain 
Tj, = Regulator stabilizing circuit time constant, 
seconds 

Sj, = Exciter saturation function 

A, 3 = Exciter saturation constants. 


The analog representation of excitation system is also shown 
in Eig.4.3, from which the dynamic equations for digital 
computer solution are written. Eor calculating initial 
conditions, the following equations are used. 


^5 - 

^3 “ 

E 


Y. 


Ref 








(4.7) 

(4.8) 

(4.9) 


The initial field voltage, E^,, is set equal to E^ calculated 
from equation (4.2). ■ IXiring step by step integration the 
following equations are used. 


^1 = \ef - ""t 
^2 " tJ ^E " ^3 


( 4 . 10 ) 


(4.11) 
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X 4 

= 

- ^2 

( 4 . 12 ) 


= ^5 

’ i ’S’ i max 

(4.13) 


= 


(4.14) 



- =^7 

( 4 . 15 ) 


II 

Hal M 
fcd ICO 

- E 

(4.16) 

^5 

II 

X - i 

=^4 

( 4 . 17 ) 

% 

^2 

• 

(4.18) 


The dot notations, X used, represents time derivatives 
of the state variables. The differential equations (4.16) 
to (4.18) are solved at each time step to give effectively by 


controlled field voltage, E^, 
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CHAPTER 5 

IMPLMEHTiiTIOF IHGLUEIHG ElKTRICAL LOAEIFG 


5 .1 General Elscussion 

In Chapter 3, we assumed that the total electrical 
demand of the loads plus network losses does not change 
even after the disturbance (generator trip-out). "We wish 
to study this assumption in greater detail. Eollowing 
generation trip-out, there will be changes in network losses 
due to redistribution of powers. Similarly there will be 
changes in voltages and frequency and this could also change 
load demand. Therefore, the total electrical power may 
change. 

¥e shall now include representations of network, 
generator electrical parameters, excitation system and 
nature of loads. If loads were represented as constant 
PjQ type, obviously ^gip would also be a constant except 
for losses in the network. We shall represent loads by 
constant admittance which is given by 



where P and are the specified real and reactive powers 
P * 

for the pth load and is the load bus voltage magnitude 
prior to the distirbance. 
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5«2 Calculation of Total Electrical Power Output 

'let 0 . be the jth rotor angle measured vith respect 
J 

to system centre of inertia. Then the rotor angle 9. measured 

J 

■with respect to synchronous reference frame can be expressed 
as 


" ®o 


(5.1) 


where 0^ = system centre of inertia or system centre of 
angle 

0. = rotor aigle of 3th machine measured with respect 
to system centre of inertia 

and 9- = rotor angle of jth machine measured with respect 
to synchronous reference frame. 

Now 3th rotor equation of motion 


M . w. = (p . - P .) 

3 D m3 63^ 

can be re-written as 

9;) = - Pej) 


(5.2) 


(5.3) 


^ *'• 

where ’9^+9^ = 9^ 


and = “n “ 


1 _ 

M, 


(1). 


(P 


^T 


mT 


^eT^ 


(5.4) 


Substituting equation (5.4) in (5.3), it becomes 


M 


= ^m3 " V^mT " ” ^3 ^3 


(5.5) 


Neglecting oscillatory power transfers from one rotor to 
another resulting from synchronizing oscillations leads 
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to the assumption that no generator rotor can he accelerating 
or decelerating relative to another during the time step of 
integration. It follows, therefore, that during this time, 
no rotor can he accelerating or decelerating relative to the 
system centre of inertia. This means that any change that 

takes place in 0^ during the time step of integration is 

« • 

assumed to occur at constant velocity. Hence 0. in 
equation (5.5) can he set equal to zero. Therefore, generator 
j must he producing electrical power as given hy the equation 
(5.6). 



M. 

M^^^T “ ^eT^ 


(5.6) 


The system equations of motion are given helow. 


d w , 

— J C p _ p ^ 

mT eT'^ 

d 9^ 

0 _ u U 'V 

dt “ ^ 0 " s^ 


(5.7) 

(5.S) 


and “o^ where j =1,2,.. .,n (5.9) 

At any instant of time, knowni'g equations (5.7) to 

(5.9) can he integrated forward hy one time step At. Prom 
this the new values of and are known at (t+ At). 

Only P m need to he calculated at (t + At) so as to take 
care of changes in network losses and load powers. At this 
stage an iterative procedure is adopted hy using a power 
flow solution method to calculate new P^g, as described 


helow. 
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Prom the initial load flow conditions obtained prior 
to disturbance, E^ and E^ are calculated using equations 
(4.1) to (4.3) respectively, Thr initial field voltage 
is set equal to E^. The excitation systan equations are 

written symbolically as 

Ej, = P(Ej„ I Vj |) (5,10) 


The equations (4.4) and (5,10) are integrated forward one 
time step for each generator and the new values of Ej 
calculated from equation (4.5) set up the necessary 
constraint over the magnitude of generator internal 
voltage E via equation (4.6) and this together with the 
generator power calculaired from equation (5.6) are specified 
at all generator internal buses which are added to the 
network, except at the slack bus where angle and are 
specified,. These buses are treated as the voltage controlled 
buses. 


The generator admittance, 1/(3^3 » becomes a 

branch of the network. The generator terminal buses are 
treated as load buses with P and Q specified as zero. The 
loads are represented as constant admittances by using the 
relation 




(5.11) 


The diagonal elements of Y-matrix are modified accordingly 
by Y . All the necessary information is available to carry 
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out a power flow solution. Gauss-Seidel power flow solution 
method was used to compute -I’gQi in the power-frequency tran- 
sient simulation program. 

Let the slack generator hus power obtained from power 
flow solution provides the power needed to 

balance total network input against loads plus losses and 
let the slack generator power calculated from equation (5.6) 
’^es' * ^slack from Lgg> it indicates that the 

slack generator is not having the same acceleration as the 
other generators which is contrary to the dynamic energy 
balance simulation procedure. However, the power flow 
solution provi4es a new value of Pg^j, and this can be expres- 
sed as 

n 

‘ ^eT(new) ” ^ejCold) ^slack (5.12) 

I 

where s = Slack machine number. 

The above equation can also be written as 


^eT(new) “ ^eT(old) ^^slack " ^es^ (5.13) 

Prom equations (5.6) and (5*13), it can be shown that 

M . 

^e3(new) “ ^e 3 (oia) slack 

where 3 = 1 ,2, . , . ,n. 

The equations (5.14) provide new generatorpowers which form 
the input to a new power flow solution. The iterative 
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procedure between equations (5.14) and power flow solution 
is continued until - ^es^ ' — ^ where e is a specified 

tolerance. At this stage, the at (t +At) is correctly 
computed. Since all the conditions at (t +a t) are known, 
the eiquations of motion of system centre of inertia can be 
integrated foiward one time step. 

5.3 Digital Simulation 

Representation of network, generators and exciters 
complicates the problem. Initially the power-frequency 
transients were simulated on IEEE-14 bus problem. A single 
line diagram of the network is given in Pig. 5.2. The bus 
conditions prior to disturbance, line data, transformer data, 
voltage controlled bus data, generators, exciters, speed 
governors and turbines data are given in ^ippendix II* 

Later on the study is carried out on UPSEB planning network 
model for which the required data is given in Appendix I. 

The data for speed governors and turbines is given in 
Tables 2.1 and 2.2, 

5 * 4 Step by Step Algorithm for Implementation on Digital 
Computer 

Step - 1 ; Read bus data corresponding to pre disturbance 

conditions and also data for lines, transformers, shunt 

loads, generators, speed governors and turbines, 

excitation systems. The generator inertia constant H, 

X - , XI and X„ are converted to the common base chosen 
d ^ d q 

(200 MVA). 
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Step - 2 : Prom the line, transformer and shunt load data, 
Y-matrix is constructed. 

Step - 3 : Generator internal buses are added to the network 
through the generator admittance 1/(r. + j X ) which becomes 

a. 

a branch to the network. These buses are treated as voltage 
controlled buses where voltage magnitude} j is specified 
by equation (4.1) and P is specified from initial bus power 
conditions. 


Step - 4: Initial conditions of state variables of 
generators, Xq, exciters, X^ and speed governors and 
turbines, X|p,are calculated from equations (4.1) to (4.3), 
(4.7) to (4. 9), and (2.1) to (2.6) (for thermal), (2.16) to 
(2.13) (for hydro) respectively. 


Step r> 5 ; Prom initial conditions Pgj* ^nd 0^ are 

calculated by making use of equations (3.1), (3.2) and (1.1) 
respectively. Total moment of inertia of the system is 
calculated by 


n 

I 

0=1 ^ 


Initially, 


1,2,..., n . 


Step - 6 ; 


Generator trip-*out condition is simulated and 


this is done as shown below.. 



= 0.0 


= 0.0 




where j=Trip-out generator 
numb er . 
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I'oad flow is carried out to obtain voltages and 
immediately after the disturbance. 

— r,. .?■ * I^ime derivatives of state variables of excitation 

« 

systems, speed governors and turbines, and generators, 
Xq. are calculated by equations (4.10) to (4.18), (2.7) to 
(2.15) (for thermal), (2.20) to (2.29) (for hydro) and (4.4) 
respectively and also 

d (Oq "I d 0 

dt' “ ^eT^ ~dt^ ^ ^ “o “ “s^ 

are formed. 


■Step - 9 ‘, The differential equations obtained from Step-8 
are solved by Modified Euler’s method of integration. At 
this stage, state variables at (t + At) are obtained. Time T 
is advanced by one time step At. Prequency f = w^/2ir is 
cal oulated. 

Step - 10 : Few values of |E | are specified at all generator 

Si 

internal buses by the equation (4.6). The generator powers, 
P, calculated from equation (5.6), are also specified at 
these buses, except at the slack bus. The slack bus angle 
is specified. 

Step - 11 : I/oad flow is performed to obtain new voltages and 
slack bus power. let it be called as . Compute 

slack machine power ’Egg’ from equation (5.6). 

Step - 12 : If (Eg^gQ^j. - Egg) < c , wh 0 r 6 e is the specified 
tolerance* is computed and enter at Step-9. Otherwise 
go to Step-15. 
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St6£_^_13: At the generator internal buses, power P. is 
modified by 

14 . 

^e3(new) = ^ejCold) ^ slack ' ^es^* 

St_e] 2 _j-__H : Steps 11 to 14 are repeated until (P , -P ) <.e . 

SXclClC GS 

At this stage all s-t (t + At) are oaloulated correctly 

and Pgj is computed. 

~ ^ 5 1 Steps S to 15 are repeated until the simulation 
period is finished. 


Stej 2 ___j-_J^: The frequency is plotted against time. Interrelation- 
ships between various subroutines is shown in ^ig.5.1. 

5«5 Computer Results for IEEE-14 Bus System ; 


General Information: 

Tfumher of buses 

# 

iNTumber of linos 

Wumber of voltage controlled buses 
Kumber of shunt loads 
N’umber of generators 
formal operating frequency 
System base MVA chosen 


= 14 

= 20 
= 4 

= 1 
= 2 

= 50.0 Hz 

= 100.0 MYA 


Total moment of inertia on 100HVA base = 0.102. 

= 514.0 rad/sec. 


Initial velocity of system centre of 

inertia ' 


Initial sys'teni centre of angle 9^ = 0.3063 radians 


Integration time step used = 0.1 sec. 

Generator at bus Ho. 2 meeting a load demand of 40.0 W out 
of 272,4 MW is tripped-out, i.e. about 15?^ generation is 
dropped. 
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The total spinning reserve available at the time of dis- 
turbance which is sufficient to overcome the generation 
loss = 52.35 M¥A. 

The total electrical power generated prior to disturbance 

= 2.724 p.u. 4m. 

The total mechanical power output prior to disturbance 

= 2.732 m. p. a. 

The simulation results for loss of generation on 
IEEE-14 bus system are given in Table 5.1, 

It is observed from the results given in Table 5.1, 
that the pre-disturbance (t = 0") electrical power is 
2.724 p.u. and immediatr-ly after the disturbance (t = 0"^) 
is 2.733 p.u. ■!#. The change in is therefore 0.34^ 

compared to a 15^ loss of generation. In 1972 Berg [ 6 1 
reported actual field tests on a Norweign system. The 
total installed generating capacity of the system exceeded 
2000 MW and amounts of generation loss ranged from 3?^ to 
of the system load. The spinning reserve was adequate 
to compensate for generation loss. The very interesting 
conclusions were that a loss of generation of only a 
few percent had very little effect on the system voltage 
and frequency. Eor as great a loss as 27^, the voltage 
drop was less than 2 percent for the HV network, indicaing 
that the total drop in electrical load is very small. This 
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is confirmed in the above IHEE-14 bus study carried out in 
time steps of 0.1 seconds upto 5.6 seconds. The reason for 
carrying out the study upto 5.6 seconds only was that in 
simulation program a not very efficient Oauss-Seidel load 
flow computer program had to be used repeatedly. It took 
26 minutes of computer time in IBM 7044 to simulate only 
5.6 seconds of transient. The graph is shown in I'ig. 5.3. 
Such a large amount of computing time may be attributed to 
possible inefficient program besides the Gauss-Seidel load 
flow program, this stage we wish to observe that it 
may be adequate to treat the electrical loading as constant 
at the pre-disturbance level. This observation is further 
substantiated when we take up the IIP State planning model 
next. 

5.6 Oomputer Results for UPSEB Planning Model Network 


General Information: 

Number of buses = 96 

Number of lines =156 

Number of voltage controlled buses = 27 

Number of shunt loads = 44 

Number of generators = 20 

Normal operating frequency = 50.0 Hz 

System base M7A chosen = 200.0 

Total moment of inertia on 200 MYA base = 0.602 


Initial velocity of systan centre of 

inertia 


= 314.0 rad/sec. 
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Initial angle of system centre of inertia = = 0.2102 rad. 

Integration time step used = 0.1 

Slack: bus number =94. 

Generator at Anpara power station having a capacity of 52S 
MYA meeting a load demand of 400 MW out of 5556.’87 MW was 
tripped-out, i.e. about 1 .2ffo generation was dropped. 

The total spinning reserve available at the time of 
disturbance = 560 WIA which is sufficient to overcome the 

loss of generation. 

The total electrical power generated prior to disturbance 
= 27.784 p.u. 

The total electrical power output immediately after the 
disturbance = 27.912 p.u. MW. 

The change in 1^^ after the disturbance = 0.1287 p.u. MW 
-.'Change in T'gj after the disturbance = 0.463^ 

Total mechanical power prior to 

disturbance = 27.841 p.u. MW 

Total mechanical power P^^ after the 

disturbance = 25.84 p.u. MW 

At 0.1 seconds after the disturbance, the system conditions 
are: 

Angle of system centre of inertia = 0.2102 rad, 

Telocity of system centre of inertia, “ ^ = 513.308 rad/sec. 

The frequency, f of the system = 49.99 Hz 

Total mechanical power output, = 25.839 p.u. MW 

Total electrical power output, P^^ = 27.912 p.u. MW 
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simulation "was carried out for only two time 
steps because of large computing time. However, it 
demonstrates our contention that electrical loading can be 
treated as constant. In that case the results obtained in 
Chapter 3 (Part - d) are adequate. 
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Tatle 5.1: Simulation Result for Loss of Generation 
on IEEE-14 Bus System,. 


s. 

?To. 

lime in 
sec . 

®o 

(rad. ) 

0 

(rad. /sec) 

(Ha) 

^mT 

(p.u.MW) 

p 

■^el 

(p.U.9ftF) 


0" 

0.3063 

314.000 

50.00 

2.7332 

2.724 


0+ 

0.3063 

314.000 

50.00 

2.3332 

2.733 

1 

0.1 

0.3063 

313.687 

49.99 

2.331 

2.733 

2 

0.2 

0.2591 

313.214 

49.92 

2.331 

2.733 

3 

0.3 

0.1646 

312.663 

49.85 

2.264 

2.733 

4 

0.4 

0.0149 

312.156 

49. ?6 

2.302 

2.733 

5 

0.5 

-0.1854 

311.583 

49.68 ■ 

2.245 

2.733 

6 

0.6 

-0.4430 

311 .014 

49.59 

2.249 

2.733 

7 

0.7 

-0.7576 

310.417 

49.59 

2.226 

2.733 

8 

0.8 

-1.1318 

309.847 

49.4-0 

2.248 

2.733 

9 

0,9 

-1 .5630 

309.243 

49. ■^I 

2.220 

2.733 

10 

1.0 

-2.0546 

308.690 

49.22 

2.263 

2.733 

11 

2.0 

-9.9095 

303.765 

48.41 

2.409 

2.733 

12 

3.0 

-21 .6670 

301,244 

47.96 

2.650 

2.733 

13 

4.0 

-34.6996 

301 .439 

47.96 

2.818 

2.733 

14 

5.0 

-46.8816 

302.659 

48.15 

% 

2.844 

2 .733 

15 

5.6 

-53.5S40 

303.453 

48,27 

2.847 

2.733 
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CHAPTER 6 

OVERALL OOHJLUSIOHS AHP SUGGEST lOHS EOR PU'TURE WORK 


6. 1 Overall Conolusions 

The hasic aim of this thesis vas to simulate power- 
frequency transients following a loss of large hlock of 
generation and also to review Stanton^ s power-frequency 
transients simulation procedure using dynamic energy 
"balance concepts. A realistic test problem had been chosen 
from one of the U.P. State electricity board planning model 
networks for the year 1984-35 to carry out the above studies. 

At first to get the 'feel^ for the problem, the 
frequency transients were simulated on the above system 
treating the entire system to have only thermal power 
generation and was represented by one equivalent machine 
of capacity equal to the total capacity of the system ani 
by one equivalent hydro machine of capacity equal to the 
total capacity of the system respectively and also assuming 
that the total electrical loading to be same for the pre and 
post disturbance conditions. The spinning reserves were 
adequate to overcome the loss of generation. Prom the 
above studies, it was concluded that the system which is 
predominantly thermal can pick up the load much faster 
than an hydro system and also it was observed that the 
maximum frequency deviation from, the normal value is much 
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less and settles to a steady state value faster than for 
an hydro system, fhis fact is because of smaller time 
constants involved in the thermal systen and also due to 
the smaller moment of inertia compared to that of hydro 
systems. 

In the next study the entire U.P.S.S.B. system was 
represented by one equivalent thermal machine with a 
capacity equal to all the thermal generation and one 
equivalent hydro machine equal to all the hydro generation. 
Since most of the thermal machines had similar time constants 
and th e same was true for most of the hydro machines, it 
could be done without too much difficulty. Again the 
electrical loading was treated as constant. As expected, 
the results were in between those for a purely thermal and 
for a purely hydro system studied earlier. Later on the 
identity of all thermal and hydro stations was retained 
but the electrical loading was treated as constant. The 
results were not too different than the case of one equiva- 
lent thermal and one equivalent hydro machine. This was 
because of the similar time constants of various machines 
as already mentioned. 

finally, the simulation was performed taking into 
account the change in electrical loading. This required 
the representation of electrical network, generators, 
excitation systems and loads (by constant admittance). 
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Because of the complexity of the problem, first lEBT' 14 bus 
syston was studied. It was found that the electrical loading 
immediately after the disturbance was only about 0.5^ 
different than the pre-fault electrical loading, and further 
the electrical loading stayed constant at the immediate 
post disturbance level through out the simulation period. 
Similar results were also obtained on UI'SEB planning model. 
Field studies carried out by Berg on the actual Norwegian 
system also show similar results. 

Thus it can be concluded that if the loads are 
treated either by constant P,Q or constant admittance 
representation, then the power-frequency transients can be 
simulated in a most economical way by considering ’ 
only speed governor and turbine actions and treating the 
electrical loading constant. Fig. 6.1 shows the system 
frequency for three cases fiz., 1) Treating entire UPSEB 
system by one equivalent thermal and one equivalent hydro 
machine, 2) Retaining the identity of all thermal and 
hydro stations and 3) The rotor angle of a generator close 
■ to the disturbance obtained from a transient stability 
program. Curves 1 and 2 are close together for reasons 
already given earlier. The trend of these two curves is 
similar to that obtained by transient stability study. 
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6.2 Suggestions for iPuture Work 

In this thesis, the. loads have heen represented by 
constant admittances. It -would be desirable to consider other 
possible models like constant current, constant MYA and 
frequency dependence etc . More detailed models for large 
induction and synchronous motors -which can take into account 
voltage and frequency changes should also be examined. 

In tie simulation method for power-f requency transients, 
the whole system was treated to have one equivalent velocity 
of system centre of inertia which may not be true. It may be 
necessary for large systems to divide into -two or three 
equivalent areas. Unless this subdivision into different 
areas is performed, the power-frequency transients which 
we obtain may be quite different from the realistic values. 



KEPERENCES 


59 


1 . K.Neil Stanton, "iDynainic Energy Balance Studies for Simulation 
of Power-Erequency Transients," Proc, Power Industry Computer 
Applications C onferenoe, Boston, Mass, May 1971. 

2. R, Podmore, "Power System I^namic Simulation Program", 

Users Manual, 

5. E.W. KimUarlc, Synchronous Machines, Power System Stability, 
’7ol.3, Hew York, Dover Publications Inc,, 196S. 

4. IEEE Committee R.eport, "Computer Representation of Excitation 
Systems", IEEE Trans. PAS, Vol. 87, pp. 1460-1464, June, 1968. 

5. Ulenn W, Stagg and Ahmed, H. El-Abiad, Computer Methods in 
Power System Analysis, Tokyo, McG-raw-Hill Kogakusha Ltd., 

1968. 

6. G.J.Berg, "Systsn and Load Behaviour Poliowing Loss of 
G-eneration, Experimental Results and Evaluation", Proc. lEE, 
Vol.119, No. 10, October 1972. 

7. V, Venikov, Transient Processes in Electrical Power Systems, 
Moscow, Mir Publishers, 1977. 

8. A, Barzam, Automation in Electrical Power Systems, Moscow, 

Mir Publishers, 1977. 

9. IEEE Press Press Publication, Stability of Large Electric 
Power Systems, New York, IEEE Press, 1974. 

10. Ollel.Elgerd, Electric Energy Systems, Theory, An Introduction, 
New Delhi, Tata McGraw-Hill Publishing Company Ltd., 1971. 

11. IEEE Committee report, "Ijynamic Models for Steam and Hydro 
Turbines in Power System Studies", IEEE Trans»,PAS, Nov/Dec., 
1973. 

12. H.E. Lokay and V. Burtnyk, "Application Underfrequency 

Relays for Automatic Load Shedding", IEEE Trans., PAS, 

January 1970, 



APPETOX I 


60 


UPSEB PLAIWITTG MODEL TO\fORK: DATA 


a) Bus Conditions Prior to Disturlpance 


Bus 

JTo. 

GENEHATIOE(MW) 

P - Q 

LOAD POVfEIl(MW) 

P Q 

Voltage Phase 

’ Magnitude angle 
in p.u. in deg. 

1 

2 

3 

4 

5 

6 

7 

1 

700.00 

155.45 

0.00 

0.00 

1 .05 

-1.528 

2 

0.00 

0.00 

0.00 

0.00 

1 .06 

-5.637 

3 

300.00 

13^^. 07 

0.00 

0.00 

1 .02 

-3.443 

4 

0.00 

0.00 

0.00 

0.00 

1.10 

-10.946 

5 

736.00 

460.00 

0.00 

0.00 

1.04 

2.934 

6 

0.00 

0.00 

0.00 

0.00 

1.07 

-5.157 

7 

63.00 

34.00 

0.00 

0.00 

1.04 

-9.123 

3 

0.00 

0.00 

40.00 

25.00 

1.09 

-12.644 

9 

197.00 

93.00 

0.00 

0.00 

1.03 

-10.117 

10 

0.00 

0.00 

230.00 

123.30 

1 .09 

-12.735 

11 

0.00 

0.00 

31.00 

53.00 

1.07 

-14.358 

12 

0.00 

0.00 

0.00 

0.00 

1.03 

-14.399 

13 

0.00 

0.00 

170.00 

126.00 

1 .06 

-17.431 

14 

0.00 

0.00 

0.00 

0.00 

0.97 

-11.163 

15 

0.00 

0.00 

0.00 

0.00 

1 .02 

-15.573 

16 

0.00 

0.00 

154.00 

134.00 

1.02. 

-20.630 

17 

0.00 

0.00 

91 .00 

67.00 

1 .00 

-21 .606 

18 

0.00 

101.42 

0.00 

0.00 

1.00 

-20.285 

19 

0.00 

0.00 

207.00 

154.00 

0.99 

-23.731 

20 

0.00 

0.00 

0.00 

0.00 

0.97 

-21.693 


C ontd. . 
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1 

2 

ry 

D 

4 

21 

0.00 

0.00 

72.00 

22 

0.00 

0.00 

52.00 

23 

0.00 

0.00 

86.00 

24 

0.00 

0.00 

0.00 

25 

0.00 

0.00 

87.00 

26 

0.00 

0.00 

0.00 

27 

0,00 

0.00 

142.20 

28 

0.00 

0.00 

77.00 

29 

190.00 

120.00 

0.00 

30 

0.00 

0.00 

49.00 

31 

0.00 

0.00 

0.00 

32 

0.00 

0.00 

235.00 

33 

95.00 

39.89 

0.00 

34 

0.00 

0.00 

27.00 

35 

0.00 

0.00 

0.00 

36 

150.00 

90.00 

0.00 

37 

0.00 

0.00 

77.00 

3 S 

0.00 

0.00 

0.00 

39 

0.00 

0.00 

227.00 

40 

0.00 

0.00 

0.00 

41 

0.00 

0.00 

0.00 

42 

0.00 

0.00 ' 

212.00 

43 

0.00 

0.00 

0.00 

44 

0.00 

0.00 

74.00 


5 

6 

7 

52.00 

1 .00 

- 24.595 

42.00 

0.96 

- 20.441 

63.00 

1 .01 

- 23.530 

0.00 

1.10 

- 11.592 

64.00 

1.06 

- 36.738 

0.00 

1.02 

- 11.546 

100.00 

1 .02 

- 17.664 

56.00 

1 .05 

- 20.671 

0.00 

0*95 

- 15.826 

35.00 

1.03 

- 17.824 

0.00 

1.02 

- 17.591 

175.00 

1 .06 

- 20.503 

0.00 

1.05 

- 25.967 

18.00 

1 .09 

- 27.757 

0.00 

1.01 

- 16.470 

0.00 

0.97 

- 14.050 

57.00 

1 .02 

- 20.235 

0.00 

0.92 

- 24.184 

170.00 

1 .02 

- 24.186 

0.00 

1.03 

- 13.219 

0.00 

1 .04 

- 18.023 

159.00 

1 .12 

^ 22.058 

0.00 

1.00 

- 25,320 

55.00 

1 .08 

- 28.925 

0 ontd . . . 
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1 

2 

3 

4 

5 

6 

7 

45 

0.00 

34.06 

0.00 

0.00 

1 .00 

- 30.522 

46 

0.00 

0.00 

73.00 

65.00 

1.05 

- 36.694 

47 

0.00 

0.00 

0.00 

0.00 

0.91 

- 41.339 

43 

0.00 

0.00 

153.00 

135.00 

0.93 

- 41.393 

49 

35 .00 

26.00 

0.00 

0.00 

1 .05 

- 37.435 

50 

0.00 

0.00 

13.00 

11.40 

1 .00 

- 41.055 

51 

176.00 

50.26 

0,00 

0.00 

1 .03 

- 31 .356 

52 

0.00 

0.00 

3.00 

5.00 

1.05 

- 36.948 

55 

0.00 

0.00 

25.00 

22.00 

1 .02 

- 39.506 

54 

0.00 ' 

2 S .43 

45.00 

39.00 

1.00 

- 42.424 

55 

0.00 

0.00 

0.00 

0.00 

1 .03 

- 33.411 

56 

0.00 

0.00 

0.00 

0.00 

0.93 

- 37.353 

57 

0.00 

0.00 

216.00 

190.00 

1.02 

- 40.325 

53 

0.00 

0,00 

139.00 

104.00 

1 .03 

- 39.599 

59 

0.00 

0.00 

241 .00 

241.00 

1 .02 

- 30.342 

60 

400.00 

243.00 

0.00 

0.00 

0.99 

- 27.560 

61 

0.00 

0,00 

0.00 

0.00 

1.04 

- 33.236 

62 

191 .00 

114.00 

0.00 

0.00 

0.96 

- 30.937 

63 

0.00 

0.00 

0.00 - 

0.00 

1.03 

- 34.291 

64 

50.00 

39.31 

0.00 

0.00 

1 .05 

- 55.448 

65 

0.00 

0.00 

151.00 

151 .00 

1.07 

- 37.811 

66 

0.00 

5’^. 35 

0.00 

0.00 

1 .00 

- 40.337 

67 

0.00 

0,00 

191.00 

169.00 

1.00 

- 44.161 


C ontd . . 
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1 

2 

3 

4 

5 

6 

7 

68 

0.00 

0.00 

0.00 

0.00 

1 .02 

-34.749 

69 

0.00 

0.00 

175.00 

153.00 

1 .05 

-38.901 

70 

0.00 

0.00 

720.00 

520.00 

0.91 

-30.837 

71 

0.00 

335.94 

32.00 

34.00 

1 .00 

-35.638 

72 

0.00 

0.00 

132.00 

132.00 

1 .01 

-39.737 

73 

0.00 

0.00 

0.00 

0.00 

0.97 

-33.304 

74 

0.00 

0.00 

172.00 

172.00 

1 .00 

-40.865 

75 

0.00 

0.00 

69.00 

69.00 

1.02 

-35.611 

76 

0.00 

0.00 

0.00 

0.00 

1 .04 

-33.422 

77 

0.00 

0.00 

102.00 

102.00 

1.07 

-36.882 

7S 

43.00 

24.00 

0.00 

0.00 

0.98 

-33.657 

79 

0.00 

0.00 

0.00 

0.00 

1.06 

-36.954 

SO 

0.00 

0.00 

0.00 

0.00 

0.99 

-36.393 

81 

0,00 

0.00 

75.00 

58.00 

1 .06 

-37.395 

82 

0.00 

0.00 

115.00 

115.00 

0.97 

-39.655 

83 

152.00 

40.09 

0.00 

0.00 

1.05 

-27.958 

34 

0.00 

0.00 

3.00 

6.00 

1 .09 

-31.759 

85 

272.00 

80.00 

0.00 

0.00 

0.93 

-25.924 

36 

0.00 

0.00 

0.00 

0.00 

1 .07 

-30.318 

37 

39.00 

45.00 

0.00 

0.00 

0.99 

-26.834 

88 

0.00 

0.00 

0.00 

0.00 

1.06 

-31 .652 

39 

191 .00 

95.00 

0.00 

0.00 

1 .00 

-26.336 

90 

0.00 

0.00 

0.00 

0.00 

1.03 

-30.449 

91 

0.00 

0.00 

0.00 

0.00 

0.96 

-33.729 


Jontd 
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12 3 4 5 6 7 

92 0.00 0.00 55.00 53.00 1 . 0 ? - 33.734 

93 0.00 0.00 62.00 46.00 1.07 - 35.530 

94 1151.37 239.33 0.00 0.00 1.03 0.00 

95 0.00 0.00 0.00 0.00 1.07 - 5.394 

96 365.00 156.19 0.00 0.00 1.03 - 12.725 
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(■b) Line Data (MVa Base = 200.0). 


Line 

ITo. 

Prom 

Bus 

To 

Bus 

Line Impedance 

Half Line 

Charg. Admit. 

(p-^») 

Off Horn. 
Tr. Turns 
Ratio 

1 

2 

3 

4 

5 

6 

7 

8 

1 

1 

2 

0.00000 

0.02380 

0.0 

0.00000 

0.950 

2 

2 

4 

0.00000 

0.06250 

0.0 

0.00000 

0.925 

3 

2 

6 

0.00123 

0.01248 

0.0 

0.04077 

1 .000 

4 

95 

14 

0.00817 

0.08320 

0.0 

0.27180 

1 .000 

5 

2 

26 

0.00972 

0.09900 

0.0 

0,32344 

1.000 

6 

3 

4 

0.00000 

0.03640 

0.0 

0.00000 

0.900 

7 

4 

S 

0.00940 

0.08303 

0.0 

0.00026 

1 ,000 

3 

4 

10 

0.06552 

0.35417 

0.0 

0.00425 

1 .000 

9 

4 

12 

0.01637 

0.07823 

0.0 

0,06669 

1 .000 

10 

4 

31 

0.01976 

0.09608 

0.0 

0.18158 

1 .000 

11 

4 

24 

0.01171 

0.05943 

0.0 

0.01180 

1 .000 

12 

5 

6 

0.00000 

0.04760 

0.0 

0.00000 

0.900 

13 

6 

35 

0.01735 

0.10506 

0.0 

0.57010 

1 .000 

14 

7 

8 

0,00000 

0.22200 

0.0 

0.00000 

0.925 

15 

S 

10 

0.03184 

0.0757S 

0.0 

0.00872 

1 .000 

16 

8 

11 

0.03184 

0.07578 

0.0 

0,00872 

1 .000 

17 

9 

10 

0.00000 

0.05700 

0.0 

0.00000 

0.925 

IS 

24 

10 

0,00000 

0.10000 

0.0 

0,00000 

1.000 

19 

10 

11 

0.05282 

0.12571 

0.0 

0.01447 

1.000 

20 

11 

13 

0.06706 

0.15959 

0.0 

, 0.01837 

1.000 

21 

12 

13 

0.00000 

0.06860 

0.0 

0.00000 

0,925 


0 ontd, . 
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1 

2 

3 

4 

5 

6 

7 

8 

22 

12 

15 

0.03514 

0.17829 

0.0 

0.03540 

1.000 

23 

13 

17 

0.09150 

0.21777 

0.0 

0.02507 

1 .000 

24 

14 

15 

0.00000 

0.06250 

0.0 

0.00000 

0.900 

25 

15 

16 

0.00000 

0.10000 

0.0 

0.00000 

0.925 

26 

15 

18 

0.03514 

0.17828 

0.0 

0.03540 

1 .000 

27 

16 

17 

0.0415S 

0.09896 

0,0 

0.01139 

1 .000 

28 

17 

19 

0.09946 

0.23671 

0.0 

0,02725 

1.000 

29 

IS 

19 

0.00000 

0.06000 

0.0 

0.00000 

0.975 

30 

18 

20 

0.02175 

0.11037 

0.0 

0.02192 

1 .000 

31 

IS 

30 

0.01332 

0.06708 

0.0 

0.05448 

1 .000 

32 

19 

21 

0.12036 

0.30234 

0.0 

0.00781 

1 .000 

33 

20 

21 

0.00000 

0.20000 

0.0 

0.00000 

0.925 

34 

21 

23 

0.1244S 

0.31345 

0,0 

0.00806 

0 

0 

0 

• 

35 

22 

23 

0.00000 

0.20000 

0.0 

0.00000 

0.925 

36 

22 

27 

0.02008 

O.IOIS'^ 

0.0 

0.02023 

1 .000 

37 

23 

28 

0.11190 

0.27929 

0.0 

0.00722 

1 .000 

38 

26 

40 

0.00597 

0.06019 

0.0 

0.19665 

1 .000 

39 

26 

27 

0.00000 

0,06250 

0.0 

0.00000 

0.950 

40 

27 

28 

0,00000 

0 . 1 0000 

0.0 

0.00000 

0.950 

41 

27 

30 

0.01165 

0.05869 

0.0 

0.04766 

1 .000 

42 

28 

32 

0.20779 

0,52192 

0.0 

0.01349 

1.000 

43 

29 

30 

0.00000 

0.03960 

0.0 

0.00000 

0.900 

44 

31 

32 

0.00000 

0.04980 

0.0 

0.00000 

0.925 

45 

31 

37 

0.03459 

0.16537 

0.0 

0.14097 

1.000 

46 

33 

34 

0.00000 

0.07920 

Oontd 

0.0 

• • 

0.00000 

0.950 
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1 

2 

3 

4 

5 

6 

7 

3 

47 

34 

59 

0.04134 

0.21225 

0.0 

0.16360 

1 .000 

43 

35 

37 

0.00000 

0.06250 

0.0 

0.00000 

0.975 

49 

35 

70 

0.00796 

0.04365 

0.0 

1 .05943 

1 .000 

50 

35 

2 

0.01592 

0.09734 

0.0 

0.53001 

1 .000 

51 

36 

37 

0.00000 

0.15340 

0.0 

0,00000 

0.900 

52 

37 

39 

0.00000 

0.06660 

0.0 

0.00000 

0.950 

53 

37 

59 

0.02594 

0.12393 

0.0 

0.10569 

1.000 

54 

33 

39 

0.00000 

0.27500 

0.0 

0.00000 

0.900 

55 

40 

41 

0.00000 

0.06250 

0.0 

0.00000 

0.950 

56 

41 

37 

0.01016 

0.04354 

0.0 

0.04133 

1 .000 

57 

41 

42 

0.00000 

0.06440 

0.0 

0.00000 

0.900 

53 

41 

43 

0.02677 

0.13534 

0.0 

0.02693 

1.000 

59 

42 

44 

0.16030 

0.40390 

0.0 

0.01044 

1.000 

60 

43 

44 

0.00000 

0.20000 

0.0 

0.00000 

0.900 

61 

43 

45 

0.02845 

0.14433 

0.0 

0.02366 

1 .000 

62 

44 

46 

0.17241 

0.43305 

0.0 

0.01119 

1.000 

63 

45 

46 

0.00000 

0.20000 

0.0 

0.00000 

0.925 

64 

46 

43 

0.07643 

0.13133 

0.0 

0.02094 

1 . 000 

65 

47 

43 

0 . 00000 

0.20000 

0.0 

0.00000 

0.925 

66 

43 

50 

0.07336 

0.13763 

0.0 

0.02160 

1.000 

67 

43 

54 

0.16343 

0.42417 

0.0 

0.01091 

1 .000 

63 

43 

57 

0.03317 

0.19793 

0.0 

0.02273 

1.000 

69 

49 

50 

0,00000 

0.33000 

0.0 

0.00000 

1.000 

70 

51 

52 

0.00000 

0.11400 

0.0 

0.00000 

1.000 


0 ontd. . . 
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1 

2 

3 

4 

5 

6 

7 

S 

71 

52 

53 

0.04680 

0.11760 

0.0 

0.01215 

1.000 

72 

52 

53 

0.05339 . 

0.12705 

0.0 

0.01462 

1 .000 

73 

53 

54 

0.09360 

0.23510 

0.0 

0.00607 

1 .000 

74 

53 

57 

0.09921 

0.24921 

0.0 

0.00604 

1.000 

75 

56 

55 

0.02295 

0.23410 

0.0 

0.06965 

1 .000 

76 

56 

61 

0.02577 

0.12964 

0.0 

0.02641 

1 .000 

77 

57 

53 

0.06359 

0.15135 

0.0 

0.01742 

1 .000 

73 

57 

72 

0 . 239 S 0 

0.60373 

0.0 

0.01553 

1 .000 

79 

53 

79 

O.OS 143 

0.20497 

0.0 

0.02108 

1 .000 

SO 

58 

SI 

0 . 0 S 148 

0.19392 

0.0 

0.02232 

1.000 

S 1 

59 

63 

0.02133 

0.10194 

0.0 

0.03690 

1.000 

32 

60 

61 

0.00000 

0.05690 

0.0 

0,00000 

0.900 

S 3 

61 

63 

0.01165 

0.05363 

0.0 

0,01194 

1 .000 

S 4 

61 

6 S 

0.00666 

0.03133 

0.0 

0.02714 

1.000 

S 5 

62 

63 

0.00000 

0.06740 

0.0 

0,00000 

0.900 

S 6 

63 

65 

0.00000 

0.07760 

0.0 

0.00000 

0.950 

S 7 

63 

66 

0.03179 

0.16131 

0.0 

0.03203 

1.000 

SS 

63 

6 S 

0.00749 

0.03531 

0.0 

0.03053 

o 

o 

o 

* 

S 9 

64 

65 

0.00000 

0.19500 

0.0 

0.00000 

0.950 

90 

65 

67 

0.13210 

0.45750 

0.0 

0 . 011 S 2 

1.000 

91 

65 

69 

0.08424 

0.21159 

0.0 

0.00547 

1 . 000 

92 

66 

67 

0.00000 

0.10000 

0.0 

0.00000 

0.975 

93 

6 S 

69 

0.00000 

0.09700 

0.0 

0.00000 

0.950 

94 

6 S 

71 

0.00999 

0.04775 

0.0 

0.04070 

1.000 
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1 

2 

3 

4 

5 

6 

7 

8 

95 

69 

72 

0.11232 

0.28212 

0.0 

0,00729 

1.000 

96 

70 

71 

0.00000 

0.06250 

0.0 

0,00000 

1.050 

97 

71 

55 

0.02436 

0.24760 

0.0 

0.07380 

1.000 

98 

71 

73 

0.01466 

0.07437 

0.0 

0.01477 

1.000 

99 

71 

72 

0.00000 

0.07800 

0.0 

0.00000 

0.975 

100 

71 

75 

0.02884 

0.14557 

0.0 

0.02942 

1 .000 

101 

72 

74 

0.02274 

0.05713 

0.0 

0.00590 

1.000 

102 

73 

74 

0.00000 

0.09320 

0.0 

0.00000 

0.950 

103 

75 

76 

0.02314 

0.11642 

0.0 

0.02371 

1.000 

104 

76 

77 

0,00000 

0.10180 

0.0 

0.00000 

0.950 

105 

76 

86 

0.01669 

0.08432 

0.0 

0.06784 

1.000 

106 

76 

90 

0.04016 

0.20376 

0.0 

0.04046 

1 .000 

107 

77 

81 

0.02875 

0.06843 

0.0 

0.00787 

1 .000 

108 

78 

79 

0.00000 

0.27500 

0.0 

0.00000 

0.900 

109 

79 

81 

0.06552 

0.16495 

0.0 

0.00424 

1,000 

110 

79 

. 25 

0.01892 

0.04612 

0.0 

0.00121 

1 .000 

111 

80 

81 

0.00000 

0.10000 

0.0 

0.00000 

0.925 

112 

80 

55 

0.01864 

0,09380 

0.0 

0.01910 

1.000 

113 

80 

82 

0.02677 

0.13584 

0.0 

0.02697 

1.000 

114 

81 

25 

0.02434 

0.06113 

0.0 

0.00632 

1.000 

115 

82 

73 

0.02677 

0.13584 

0.0 

0.02697 

1 .000 

116 

83 

84 

0.00000 

0.10480 

0.0 

0.00000 

0.950 

117 

84 

86 

0.02434 

0.26126 

0.0 

0.00157 

1 .000 

118 

84 

92 

0.01965 

0.04997 

0.0 

0.02036 

1 .000 


0 ontd. > 
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1 

2 

3 

4 

5 

6 

7 

3 

119 

S5 

36 

0.00000 

0.06580 

0.0 

0.. 00000 

0.900 

1 20 

36 

55 

0.02664 

0.13410 

0.0 

0.02727 

1 .000 

121 

S7 

33 

0.00000 

0.21930 

0.0 

0.00000 

0,900 

122 

33 

55 

0.01673 

0,08490 

0.0 

0.06744 

1,000 

123 

39 

90 

0.. 00000 

0.03900 

0.0 

0.. 00000 

0.900 

124 

90 

55 

0.02003 

0.10133 

0.0 

0.. 0309 2 

1.000 

125 

91 

92 

0.00000 

0.20000 

0.0 

0-. 00000 

0.900 

126 

92 

93 

0.02465 

0.05960 

0.0 

0.01332 

1.000 

127 

55 

93 

0.00000 

0.10160 

0.0 

0 . 00000 

0.925 

123 

93 

25 

0.01435 

0,03613 

0..0 

0.00336 

1 .000 

129 

56 

57 

0.00000 

0,06350 

0.0 

0.00000 

0.925 

130 

96 

41 

0.00000 

0.05690 

0..0 

0,00000 

0.950 

131 

94 

95 

0.00000 

0.01390 

0,0 

0.00000 

0.950 

132 

95 

6 

0.00031 

0.00332 

0.0 

0.02700 

1.000 

133 

95 

2 

0.00031 

0.00332 

0.0 

0.02700 

1.000 

134 

95 

26 

0.00972 

0.09900 

0.0 

0.32344 

1 .000 

135 

^ 95 

35 

0.01726 

0.10500 

0.0 

0.57010 

1.000 

136 

95 

40 

0,01365 

0.14390 

0.0 

0.49000 

1.000 


71 


c) Voltage Oontrolled Bus Data 


S.No. 

Bus Ko. 

Q-Minimum 

MVAR 

Q-Maxifflum 

MVAR 

Scheduled 

hus voltageCp.u. ) 

1 

2 

3 

4 

5 

1 

3 

0.0 

214.0 

1.02 

2 

5 

0.0 

460.0 

1.05 

3 

7 

0.0 

34.0 

1.05 

4 

9 

0.0 

98.0 

1.05 

5 

29 

0.0 

120.0 

1.05 

6 

33 

0.0 

57.0 

1 .05 

7 

36 

0.0 

90.0 

1.05 

S 

38 

0.0 

100.0 

1 .05 

3 

49 

0.0 

26.0 

1.075 

10 

51 

0.0 

33.0 

1 .075 

11 

60 

0.0 

243.0 

1 .05 

12 

,62 

0.0 

114.0 

1 .05 

15 

64 

0,0 

40.0 

1.05 

14 

73 

0.0 

24.0 

1.05 

15 

33 

0.0 

76.0 ■ 

1 .05 

16 

85 

0.0 

80.0 

1.05 

17 

37 

0.0 

45.0 

1 .05 

13 

89 

0^0 

95.0 

1.05 

19 

13 

0,0 

500.0 

1.0 

20 

45 

' 6.0 

500.0 

1.0 


0 ontd. , 
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d) Shunt I^oad Data 


S.T!To. 

Bus No. 

Shunt Load AdmittanceCp.u. ) 

1 

2 

-3 


1 

11 

0.0 

0.065 

2 

13 

0.0 

0.020 

3 

16 

0.0 

0.010 

4 

17 

0.0 

0.085 

5 

19 

0.0 

0.160 

6 

21 

0.0 

0,035 

7 

23 

0.0 

0.080 

8 

27 

0.0 

0.120 

9 

28 

0.0 

0.070 

10 

32 

0.0 

0.030 

11 

39 

0.0 

0.070 

12 

42 

0.0 

0.125 

13 

44 

0.0 

0.100 

14 

46 

0.0 

0.180 

15 

48 

0.0 

0.350 

16 

54 

0 

• 

0 

0,080 

17 

57 

0.0 

0.500 

18 

58 

0 

• 

0 

0.230 

19 

59 

0 

• 

0 

0.800 

20 

65 

0.0 

0.450 

21 

67 

0.0 

0.520 

22 

69 . 

0.0 

0.450 


0 ontd 
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1 

2 

3 

4 

23 

71 

0.0 

0.250 

24 

72 

0.0 

0.500 

25 

74 

o 

• 

O 

0.520 

26 

75 

0.0 

0.280 

27 

77 

0.0 

0.270 

28 

81 

0.0 

0.120 

29 

82 

0.0 

0.520 

30 

92 

0.0 

0.030 

31 

93 

0.0 

0.020 

32 

25 

0.0 

0.050 

33 

2 

o 

9 

o 

- 0.750 

34 

6 

0.0 

- 0.250 

35 

14 

0.0 

- 0.250 

36 

26 

0.0 

- 0.500 

37 

35 

0.0 

- 1.000 

3 S 

40 

0.0 

0.000 

39 

70 

0.0 

- 0.500 

40 

95 

0.0 

- 0.500 

41 

18 

0.0 

0.500 

42 

45 

0.0 

0.420 

43 

66 

0.0 

0.265 

44 

70 

0.0 

1 .670 
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e) Generator Parameters 


S.Ho. 

Base 

MVA 

H 

(seconds) 

^a 



I 

9 

1 

940.0 

2.20 

0.00179 

2.117 

0.287 

2.117 

2 

422.5 

2.71 

0.00138 

1 .340 

0.242 

1.780 

*3 

940.0 

2.20 

0.00179 

0.287 

0.2S7 

0.287 

4 

110.0 

3.50 

0.00600 

0.901 

0.290 

0.559 

5 

335.0 

3.99 

0.00500 

0.950 

0.335 

0.600 , 

6 

275.0 

2.00 

0.00100 

2.170 

0.300 

2.050 

7 

139.0 

2.00 

0.00100 

2.170 

0.300 

2.050 

S 

217.5 

2.20 

0.00270 

1.700 

0.290 

1 .650 

* 9 

51.75 

4.00 

0.00500 

0.230 

0.230 

0.280 

*10 

220.0 

4.94 

0.00670 

0.250 

0.250 

0.250 

*11 

528. 0 

1 .91 

0.00100 

0.230 

0.230 

0.230 

*12 

340.0 

2.53 

0.00270 

0.212 

0.212 

0.212 

*13 

75.0 

3.50 

0.00500 

0.250 

0.250 

0.250 

14 

160.0 

4.55 

0.00300 

0.950 

0.300 

0.600 

15 

211 .0 

3.40 

0.00500 

0.950 

0.310 

0.600 

16 

376.0 

3.79 

0.00314 

0.950 

0.280 

0.600 

17 

103.5 

2.85 

0.00279 

0.953 

0.263 

0.600 

18 

223.0 

3.00 

0.00300 

0.950 

0.300 

0.600 

+19 

1420.0 

2.10 

0.00230 

a. 350 

0.280 

2.150 

20 

490.0 

2.20 

0.00230 

2.120 

0.287 

2.120 


*These are treated as constant excitation machines. 
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f ) Excitation System Ifata 


s. 

No . 

Gen. 
¥o • 

-A 

^E 

% 


®E 


^Rmin 

V 

Bmax 

1 

1 

250 

1.0 

0.12 

0.04 

O 

* 

o 

^.5 

-5 

5 

2 

2 

200 

1 .0 

0.175 

0.04 

1 .0 

3.5 

-5 

5 

3 

4 

200 

1 .0 

0.12 

0.50 

0.7 

2.0 

-4 

4 

4 

5 

200 

1 .0 

0.12 

0.50 

0.7 

2.0 

-4 

4 

5 

6 

200 

1.0 

0.10 

0.02 

1.33 

3.0 

-5 

5 

6 

7 

200 

1.0 

0.10 

0.02 

1 .38 

3.0 

-5 

5 

7 

3 

200 

1.0 

0.10 

0.02 

1 .38 

3.0 

-5 

5 

3 

14 

200 

1.0 

0.12 

0.50 

0.70 

2.0 

-4 

4 

9 

15 

200 

1.0 

0.12 

0.50 

0 70 

2.0 

-4 

4 

10 

16 

200 

1.0 

0.12 

0.50 

0.70 

2.0 

-4 

4 

11 

17 

200 

1 .0 

0.12 

0.50 

0.70 

2.0 

-4 

4 

12 

18 

200 

1 .0 

0.12 

0.50 

0.70 

2.0 

-4 

4 

13 

19 

600 

1 *0 

0.00 

10.72 

0.47 

100.0 

2.1 

11.5 

14 

20 

250 

1 .0 

0.12 

0.04 

0.20 

3.5 

-5 

5 
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APPENDIX II 
IEEE-14 BUS SYSTEM DATA 
a) Bus Conditions Prior to Disturbance 


Bus 

No. 

Gen era ti on (l/nv ) 
P Q 

load 

P 

Power (ITW) 

Q 

Voltage 
Magnitude 
in p.u. 

Phase 
angle in 
degrees 

1 

232.4 

-16,9 

0.0 

0.0 

1 .06 

0.0 

2 

40.0 

42.4 

21.7 

12.7 

1 .05 

-4.93 

3 

0.0 

23.4 

94,2 

19.0 

1 .01 

- 12.72 

4 

0.0 

0,0 

47. B 

-3.9 

1 .02 

-10.33 

5 

0.0 

0.0 

7.6 

1.6 

1 .02 

-3.78 

6 

0.0 

12.2 

11.2 

7.5 

1 .07 

-14.22 

7 

0.0 

0.0 

0.0 

0.0 

1.06 

-13.37 

S 

0.0 

17. S 

0.0 

0.0 

1.09 

-13.36 

9 

0.0 

0.0 

29.5 

16.6 

1.06 

-14.94 

10 

0.0 

0.0 

9.0 

5.8 

1.08 

-15.14 

11 

0.0 

0.0 

3.5 

1.3 

1 .06 

-14.79 

12 

0.0 

0.0 

6.1 

1 .6 

1 .06 

-15.07 

13 

0.0 

0.0 

13.5 

5.8 

1.05 

-15.16 

14 

0,0 

0.0 

14.9 

5.0 

1 .04 

-16.04 



"b ) Line Data 

(Base 

MYA = 100 

.0) 
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Line 

No. 

Prom 

Bus 

To 

Bus 

Line Impedance 
(p.u. ) 

Half Line 
GLarg. Admit, 
(p.u.) 

Off Norn. 

Tr .Turns 
Ratio 

1 

1 

2 

0.01938 

0.05917 

0.0 

0.0264 

1.000 

2 

1 

5 

0.05403 

0.22304 

0.0 

0.0264 

1 .000 

5 

2 

3 

0.04699 

0.19797 

0.0 

0.0219 

1 -000 

4 

2 

4 

0.05311 

0.17632 

0.0 

0.1379 

1.000 

5 

2 

5 

0.05695 

0.17338 

0.0 

0.0170 

1 .000 

6 

3 

4 

0.06701 

0.17103 

0.0 

0.1730 

1 .000 

7 

4 

5 

0.01335 

0.04211 

0.0 

0.0064 

1.000 

8 

4 

7 

0.00000 

0.02091 

0.0 

0.0000 

0.978 

9 

4 

9 

0.00000 

0.55613 

0.0 

0.0000 

0.969 

10 

5 

6 

0.00000 

0.25202 

0.0 

0.0000 

0.932 

11 

6 

11 

0.09480 

0.19392 

0.0 

0.0000 

1.000 

12 

6 

12 

0.12291 

0.25581 

0.0 

0.0000 

1.000 

13 

6 

13 

0.0661 5 

0.13027 

0.0 

0.0000 

1 .000 

14 

7 

8 

0.00000 

0.17615 

0.0 

0.0000 

1 .000 

15 

7 

9 

0.00000 

0.11001 

0.0 

0.0000 

1 .000 

16 

9 

10 

0.03131 

0.08450 

0.0 

0.0000 

1 .000 

17 

9 

14 

0.12711 

0.27033 

0.0 

0.0000 

1 .000 

18 

10 

11 

0.03205 

0.19207 

■ 0.0 

0.0000 

1 .000 

19 

12 

13 

0,22092 

0.19983 

0.0 

0.0000 

1.000 

20 

13 

14 

0.17093 

0.34302 

0.0 

0.0000 

1.000 
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c ) Voltage Gontrolled Bus Ifata 


S.TTo. 

Bus, '^ 0 , 

Q-Minimum 

MVAR 

Q-Maxiiiium 

MVAR 

Scheduled 
bus voltage 
(p.u.) 

1 

2 

-40-0 

50.0 

1 .045 

2 

3 

0.0 

40.0 

1 .010 

3 

6 

-6.0 

24.0 

1.070 

4 

8 

-6.0 

24.0 

1 .090 


d ) Shunt 

load Bata 



S.Ro. 

Bus No. 

Shunt load 

Admittance (p.u.) 

1 

9 

0.0 

0.19 


e) generator Parameters 


Base 

MVA ( ! 

H 

seconds) 

^a 

xXd 




1 

335.0 

3.99 

0.005 

0.95 


0.355 

0.60 

2 

75.0 

3.50 

0.005 

0.25 


0.250 

0.25 

f) 

Excitation 

System Data 





s. 

1^0 

Gen . TT 

. ^0. 

% 




^Rmin 

Y 

Rmax 

1 

1 200 

1.0 

0.12 0.5 

0.7 

2.0 

-4 

4 

2 

2 200 

1.0 

0.12 0.5 

0.7 

2.0 

-4 

4 



g) Turbine Parameters 
Machine No. 1 : 


80 


p 

base 

= 335.0 


= 1.0 


= 0.015 


= 5.0 


= 0.04 


= 0.30 

• 

®niin 

= -0.2 

« 

^max 

= 0.1 

%in 

= 0.0 

^ax 

= 0.B5 

Machine No. 

2; 

■p 

base 

= 75.0 

'^HP 

= 0.2 


= 4.0 

^LP 

= 1000.0 

^HP 

= 0.62 

Ppp 

= 0.38 

^IP 

= 0.0 

R 

= 0.04 


= 0.425 

* 

^in 

= -1.18 

• 

^ax 

= 0.59 

^in 

= 0.0 

®!iiax 

= 0.74 



